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Zinc and cadmium clusters interacting with a Gupta potential have previously been identified as 
prototypical metallic systems that exhibiting disordered cluster structures. Here, putative global 
minima of the potential energy have been located for these clusters for all sizes up to N < 125. 
Although none of the usual structural forms are lowest in energy and many of the clusters have no 
overall order, strong structural preferences have been identified. Many of the clusters are based on 
distorted oblate Marks decahedra, where the distortion involves the bringing together of atoms on 
either side of a re-entrant groove of the Marks decahedron. 



I. INTRODUCTION 

There has been much recent theoretical interest in the 
possibility that small clusters could have lowest-energy 
structures that are disordered or amorphous. Examples 
have been found for model clusters interacting with long- 
ranged pair potentials? 1 ^ and a wide variety of metal 
clustersi 3 i 4 i 5 i 6 i 7 i 8 i 9 i 10 . n i 12 i 13 i 14 i 15 i 16 For some of these ex- 
amples, these disordered clusters appear at sizes be- 
tween the magic numbers for the usual icosahedral, dec- 
ahedral or face-centred-cubic (fee) forms, that are typi- 
cally most stable for materials that are close-packed for 
bulk^iSii 3 ^ 1 ^* 1 ^ i.e. the disordered structures are more sta- 
ble than ordered structures with incomplete outer layers. 
More interesting are the more extreme examples, where 
disordered structures are lowest in energy even at the 
magic numbers of the common structural forms^* 7 ^ 
These results would suggest for these clusters that the 
disordered structures are dominant for the relevant size 
ranges. 

By disordered it is usually meant that the structure 
has no discernible overall order. Of course, there is still 
local order, as evidenced by structural probes, such as the 
radial distribution function, but these usually have forms 
similar to that for bulk liquids and glasses. Additionally, 
another common features for these clusters is that there 
are many other disordered structures that lie very close 
in energy to the global minimum. 

Disordered clusters are most likely to occur when the 
constraints on the nearest-neighbour distances are weak. 
Then, the energetic cost for the strain present in the 
individual nearest-neighbour bonds can be low enough 
to be offset by other advantageous features of the dis- 
ordered clusters, such as a low surface energy. For the 
Morse clusters, this occurs when the potential is long- 
ranged and has a wide, soft welli^ For metal clusters, 
the many-body part of the potential is relatively insen- 
sitive to disorder in the nearest-neighbour distances^ If 
the contraction at the surface of the cluster (again due to 
the many-body forces) is strong, the ordered structures 
can be sufficiently destabilized to make the disordered 
structures lowest in energy^ 

For most of the examples where disordered structures 
have been found to be lower in energy than the most sta- 



ble icosahedral, decahedral and fee clusters, only a few 
sizes have usually been considered So, the possibil- 
ity remains that at some other sizes the local structural 
preferences present in the disordered clusters could be 
assembled in a way that gives rise to a particularly sta- 
ble structure with overall order. If this were the case, 
global optimization over a complete range of sizes would 
then reveal these new structural motifs and magic num- 
bers. This has recently been done for model lead clusters, 
where, although none of the usual forms are ever lowest 
in energy and many of the clusters have no apparent over- 
all order, particularly stable high-symmetry clusters were 
found at some sizes iAi 

Indeed there are general grounds for expecting high 
symmetry structures to be more prevalent among the 
global minima, because such structures are likely to have 
more extreme values of the energy (both high and low). 18 
This expectation, in my experience at least, seems to be 
born out empirically. It is rare that the global optimiza- 
tion of clusters does not reveal ordered high symmetry 
forms at some sizes. For example, even for potentials 
which have been designed to favour glassy configurations, 
global optimization has revealed the presence of unusual, 
but nevertheless ordered, structuresii^i^S 

Here, I wish to examine some prototypical metal clus- 
ters that so far have only been found to exhibit disordered 
structures. For potentials of the Gupta form the depen- 
dence of the tendency to disorder on some of the param- 
eters of the potential has been elaborated^ 1 " 7 * Of the pa- 
rameterized metals zinc and cadmium clusters emerged 
as those with the strongest preference for disorder. This 
tendency was born out in global optimization studies at 
sizes where particularly stable fee, decahedral and icosa- 
hedral structures were possible; as expected, the resulting 
structures appeared disordered^ 

In this paper, I will attempt to identify structural pat- 
terns for these two cluster systems, in particular search- 
ing for new magic numbers and novel types of order. This 
aim is in a similar spirit to Ref. l2ll where possible struc- 
tural patterns for gold clusters modelled by the Gupta 
potential have been suggested. To achieve this I have 
performed global optimization for all clusters with up to 
125 atoms. A further aim is to then relate back the iden- 
tified structural patterns to the form of the potential. 
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TABLE I: Parameters for the Gupta potentials of Zn and Cd. 





P 


q 


A j eV 


e/eV 


Zn 


9.689 


4.602 


0.1477 


0.8900 


Cd 


10.612 


5.206 


0.1420 


0.8117 



II. METHODS 



A. Potential 



To model the zinc and cadmium clusters I use a Gupta 
potential^ fitted by Cleri and RosatoiSi The potential 
energy is given by 



F -^pair "i" ^cmbed 



= E<h^)+E f (^)> 



(i) 

(2) 



l<3 



where (f>(r) is a short-ranged pair potential, U{p) is a 
many-body embedding (or glue) function and pi is de- 
fined as 



P; 



= E^( r y)' 



(3) 



where p{r) is an "atomic density" function. 
For potentials of the Gupta form 

<j>{r) = 2Ae- p{r/ra ~ 1) 
F(p) = 
p(r) 



_ e -2g(r/r -l)_ 



(4) 
(5) 
(0) 



These forms arise from the second moment approxima- 
tion of a tight-binding Hamiltonian. However, these func- 
tions are non-unique. Functions that give exactly the 
same energy can be constructed by the transformation 



<t>' (r) 
F'(p) 



(f> (r) + 2gp(r) 
F(p)-9P- 



(7) 

(8) 



This transformation redistributes the total energy be- 
tween impair and -Eembcd- When 



dF 
dp 



(9) 



F'(p) has a minimum at p x tai, where p x tai is the value 
of p in the equilibrium crystal. This choice is called the 
effective pair format, and has been suggested as the most 
natural way to partition the energy between the pair and 
many-body contributions In this format, when p = 
Pxtai the pair potential controls the energy change for any 
change of configuration that does not significantly alter 
p, hence the name. More specifically, by performing a 
Taylor expansion about this reference density, one can 
show that to first order, the change in energy is due to 



the pair potential. Consequently, it is also a much more 
helpful format for relating the structure to the form of 
the potential. 

The Gupta potential in this effective pair format be- 
comes 



eff (r) = 2Ae- p{r/n '- 1) 



-2q{r/r -l) 



F eS (p) = 1 



Vp- 
p 

2 V Pxtal 



xtal 



1 



(10) 



11) 



'/'off {f) is a sum of two exponentials, so for p > 2q it is 
repulsive at short range and has an attractive well. The 
minimum in <p e s(f) is at 



r 1 + 



_J_ (4pVp** 

p-2q H 



and is of depth 



-A 



eg 



\Apyfp 



'xtal 



p/p-2q 



p~2q 



(12) 



(13) 



As p — > 2q from above, the depth of the well in the pair 
potential goes to zero. F c s(p) is quadratic in y/p and has 
a minimum of depth -£VPxtai/2 at p = p xta i- 

The Gupta parameters for the Zn and Cd potentials 
are given in Table [IJ and the functions e g, p and F e s 
are shown in Fig. The shallowness of the attractive 
well in the effective pair potential of cadmium is particu- 
larly apparent — it is only 0.26% of the depth of the min- 
imum in the embedding function. The pair well depth 
is somewhat larger for zinc, but it is still only 0.70% of 
the embedding function minimum. For both systems the 
shallowness arises because p/q is close to 2. 

This feature of the potentials will have similar struc- 
tural consequences for both systems. Firstly, the major- 
ity of the binding energy will come from the many-body 
interactions. Secondly, the pair potential provides rela- 
tively little constraint on the pair distances, except for 
the repulsion at short distances. Therefore, the most 
important feature for obtaining a low-energy configura- 
tion is to have the individual pi values as close as pos- 
sible to the optimal value. By contrast, there is little 
energetic advantage in having the nearest-neighbour dis- 
tances close to the minimum of the pair potential. For a 
cluster the former can sometimes be more easily achieved 
through a disordered structure than one based on a lat- 
tice, because the additional flexibility of not having well- 
defined nearest-neighbour distances makes it easier to ob- 
tain close to optimal p values for the surface atoms. This 
is the source of the tendency to disorder for zinc and 
cadmium clusters described by the current potentials. 

For bulk both zinc and cadmium are hexagonal close- 
packed with p^ al = 8.042, and p^ = 8 638 _ These 
are anomalously low values for close-packed materials. 
For example, for an ideal close-packed crystal with all 
nearest- neighbours at ro, the nearest- neighbour contri- 
bution to p is 12. However, the ratios of the unit cell 
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FIG. 1: The three functions that make up the potential: (a) 
<j> e a(r), p(r) and (b) F e a(p). 



parameters c/a for the Zn and Cd crystals are particu- 
larly large, and so the contribution to p from nearest- 
neighbour distances with a component in the c direction 
is significantly reduced. 

The low values of /5 x tai have structural consequences for 
the clusters. An atom can achieve the optimal p values 
with only eight or nine neighbours, so this is possible 
for surface atoms. By contrast, the nearest-neighbour 
distances for atoms in the interior need to be elongated to 
prevent unfavourably large values of p. This is bad news 
for cluster structures, such as the Mackay icosahedron, 
where the interior distances are naturally shorter than 
those on the surface. 

The exponential nature of p{r) means that significant 
changes to the p values can be achieved by relatively 
small changes to the nearest-neighbour distances. Con- 
sequently, it is somewhat easier for atoms, even those 
with low coordination numbers, to obtain nearly opti- 
mal p values. This behaviour contrasts with other many- 
body potential, for instance those produced by the force- 
matching method^ where there is no presumed form for 
p(r), and so more long-ranged functions can result 



The current potentials were obtained by keeping the 
c/a value fixed at the experimental value. This is at some 
cost in terms of the quality of the fit to other quantities. 
This led Cleri and Rosato to also construct potentials 
in which the c/a ratio was allowed to vary in the fitting 
procedure, and resulted in a better quality fit for many 
properties^ Such a potential is available for cadmium 
and has dramatically different structural properties be- 
cause p/q — 3.49. 



B. Global Optimization 

The global optimization of the zinc and cadmium 
clusters was performed using the basin-hoppin g 28 ! 29 (or 
Monte Carlo minimization^) approach. This method 
has proved particularly successful in locating putative 
global minima for a wide variety of cluster systems. The 
nature of the potentials for the current clusters, partic- 
ularly that p — 2q is small, makes global optimization 
difficult compared to many other metal clusters of sim- 
ilar size. Therefore, a considerable computational effort 
was required to extend the results up to 125 atoms. It 
proved particularly important to supplement the stan- 
dard unbiased runs from random starting points for each 
size, with short runs started from low-energy minima of 
nearby sizes with the appropriate number of atoms added 
or removed. As the structures of the two clusters are 
similar — seventy of the global minima are the same — it 
also proved useful to reoptimize the low-energy minima 
found for one metal for the other. These two types of 
'seeded' runs were applied iteratively until no further new 
global minima were located. 

It should be noted that, of course, there is no guaran- 
tee that I have been able to locate the true global min- 
ima, and the probability that a global minimum has been 
missed will increase with cluster size, as the size of the 
search space, and hence the number of minima^i^SiS 
increases exponentially with N. Examination of the 
statistics of how often independent runs locate the same 
lowest-energy minimum and the importance of the seeded 
runs can provide one with an idea of the likelihood that 
the true global minimum has been found. For example, 
for zinc clusters with less than ninety atoms virtually 
all the putative global minima were located in unbiased 
runs, but for N > 100 the majority of putative global 
minima were only located in seeded runs. Similarly, for 
the cadmium clusters virtually all the putative global 
minima were only located in seeded runs for N > 65. 
This greater difficulty is because 2p — q is closer to zero 
for cadmium. However, given the similarity of the ob- 
served structures for the two metals, I am confident that 
the vast majority of the putative global minima up to 
N = 100 cannot be bettered, but beyond this size one's 
degree of scepticism about the success of the global opti- 
mization should increase rapidly. 
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TABLE II: Energies (in eV) and point groups of putative Znjv global minima 



TV 


PG 


Energy 


TV 


PG 


Energy 


TV 


PG 


Energy 


TV 


PG 


Energy 


TV 


PG 


Energy 


3 


D 3h 


-3.751296 


28 


D 6d 


-37.076416 


53 


Cx 


-70.622509 


78 


Cx 


-104.238674 


103 


Cx 


-137.902829 


4 


T d 


-5.099346 


29 


c 2 


-38.410947 


54 


Cx 


-71.969642 


79 


c 2v 


-105.595357 


104 


Cx 


-139.255744 


5 


D sh 


-6.422287 


30 


Ci 


-39.768212 


55 


Cx 


-73.310032 


80 


c 2v 


-106.927940 


105 


Cx 


-140.593303 


6 


o h 


-7.765873 


31 




-41.107234 


56 


Cx 


-74.649595 


81 


C s 


-108.284296 


106 


Cx 


-141.928489 


7 




-9.088519 


32 


d 


-42.444864 


57 


Cx 


-75.992628 


82 


C s 


-109.625208 


107 


Cx 


-143.279843 


8 


D 2d 


-10.409083 


33 


Cx 


-43.776617 


58 


Cx 


-77.336947 


83 


C s 


-110.976689 


108 


Cx 


-144.635048 


9 


C 2v 


-11.733002 


34 


c 3 


-45.125218 


59 


Cx 


-78.672227 


84 


C s 


-112.322650 


109 


Cx 


-145.978993 


10 




-13.055156 


35 


c s 


-46.472931 


60 


Cx 


-80.019432 


85 


Cx 


-113.668080 


110 


Cx 


-147.337545 


11 


c 2 


-14.379679 


36 


D 2h 


-47.814946 


61 


c 2 


-81.363067 


86 


c 2v 


-115.021815 


111 


c 2 


-148.697585 


12 


c 2 


-15.706971 


37 


Cx 


-49.144718 


62 


Cx 


-82.708058 


87 


c 2v 


-116.375101 


112 


Cx 


-150.033337 


13 


c s 


-17.032142 


38 


c s 


-50.499983 


63 


Cx 


-84.061400 


88 


c 2v 


-117.726908 


113 


Cx 


-151.385243 


14 


Cg v 


-18.369312 


39 


c 2v 


-51.846864 


64 


c 2 


-85.419385 


89 


C s 


-119.057698 


114 


Cx 


-152.720206 


15 


c s 


-19.697479 


40 


c s 


-53.182148 


65 


Cx 


-86.749702 


90 


C s 


-120.405576 


115 


c 2 


-154.071651 


16 


c s 


-21.042019 


41 


c s 


-54.521834 


66 


Cx 


-88.097872 


91 


c 2v 


-121.767147 


116 


Cx 


-155.420956 


17 


C 2v 


-22.395676 


42 


D 4 


-55.874358 


67 


Cx 


-89.426800 


92 


c 2v 


-123.118340 


117 


Cx 


-156.762711 


18 


Cav 


-23.735972 


43 


C\ 


-57.213737 


68 


Cx 


-90.782264 


93 


C s 


-124.451825 


118 


Cx 


-158.109816 


19 


D 4d 


-25.067563 


44 


c s 


-58.554345 


69 


Cx 


-92.135925 


94 


C s 


-125.796431 


119 


Cx 


-159.451202 


20 


c 2v 


-26.395373 


45 


c 2v 


-59.895318 


70 


Cx 


-93.473840 


95 


C s 


-127.142167 


120 


c 2v 


-160.812498 


21 


c s 


-27.715028 


46 


c s 


-61.227348 


71 


Cs 


-94.828826 


96 


c 2v 


-128.501419 


121 


Cx 


-162.152965 


22 


c 2v 


-29.042499 


47 


C s 


-62.580198 


72 


C s 


-96.168254 


97 


C s 


-129.836679 


122 


Cx 


-163.502185 


23 


Cx 


-30.386460 


48 


c s 


-63.916581 


73 


Cs 


-97.513085 


98 


c 2v 


-131.170669 


123 


c 2v 


-164.859277 


24 


Co 


-31.728122 


49 


c 2v 


-65.269233 


74 


Cs 


-98.853247 


99 


C s 


-132.525475 


124 


Cx 


-166.196591 


25 


c 2 


-33.075590 


50 


C s 


-66.601903 


75 


c 2v 


-100.203580 


100 


c 2v 


-133.883827 


125 


Cx 


-167.548063 


26 


c 2 


-34.412179 


51 


Cx 


-67.938575 


76 


Cx 


-101.545341 


101 


C s 


-135.219057 








27 


C 2 


-35.744950 


52 


c 2v 


-69.281361 


77 


C s 


-102.900278 


102 


Cx 


-136.566429 









III. GLOBAL MINIMA 

The energies and point groups for the putative global 
minima are given in Tables [H] and IIIII Point hies are 
available online at the Cambridge Cluster Database. 34 
The energies of the global minima are represented in Fig. 
121 in such a way that makes particularly stable clusters 
stand out. 

The previous results for these clusters were for a se- 
lection of sizes that often show highly symmetric struc- 
tures, namely TV = 13, 38, 55 and 147. 7 For TV=13 
and 38 the three lowest-energy minima reported in Ref. 
agree with the current results. However, for Zn^, 
Zn75 and Cd75 the lowest-energy structures reported by 
Michaelian et al. correspond at best to the third, eighth, 
and twenty-second lowest-energy minima, respectively, 
and lie 0.0031, 0.0186 and 0.0086 eV above the lowest- 
energy minima reported here. These energies are signifi- 
cant compared to the variations shown in Fig. [21 I did not 
systematically attempt to optimize 147-atom clusters, as 
locating the true global minimum for this size would be 
extremely difficult. However, short basin-hopping runs 
did find structures that were 0.0966 and 0.0304 eV lower 
in energy for Zni47 and Cdi47, respectively, than the 



lowest-energy structures found by Michaelian et al.' 

Before, we examine the observed structures, there are 
a number of interesting features evident from Fig. [21 The 
energy zero in these figures is E avci a four-parameter fit 
to the energies of the global minima, where the first two 
terms correspond to volume and surface energies. For 
these two clusters the surface term is exceptionally small. 
The ratio of the surface to the volume coefficient is 7.8% 
for Zn and 3.3% for Cd. For comparison the value of 
this ratio is 197% for Lennard- Jones clusters^ 48% for 
Gupta lead clusters f 3 ^** 3 ^ 60% for aluminium clusters^ 
and 93%, 87% and 63% for Sutton-Chen silver, nickel 
and gold clusters, respectively 10 . One expects the sur- 
face energy to be lower for metal clusters than for a clus- 
ter interacting with a pair potential, because the lower 
coordinate surface atoms of a metal can increase their 
many-body embedding energy by shortening their near- 
est neighbour distances. But the Zn and Cd surface en- 
ergies are very low even for metals, because the optimal 
value of p is small enough to be achievable by the surface 
atoms and the effective pair potential only makes a small 
contribution to the energy. 

The second interesting feature is the magnitude of the 
fluctuations about the average energy. For Cd and Zn 
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TABLE III: Energies (in eV) and point groups of putative Cdjv global minima. Those labelled with a star have the same 
structure as the Znjv global minimum. 



N 


PG 


Energy 




N 


PG 


Energy 




TV 


PG 


Energy 




N 


PG 


Energy 




N 


PG 


Energy 




3 


D 3h 


-3.391649 




28 


Ci 


-32.354453 




53 


Ci 


-61.383718 




78 


Ci 


-90.441924 


* 


103 


Ci 


-119.516811 


* 


4 


T d 


-4.556244 




29 


Ci 


-33.511796 




54 


Ci 


-62.547483 




79 


Civ 


-91.609644 


* 


104 


Ci 


-120.682922 


* 


5 


D 3h 


-5.710992 


* 


30 


c 2 


-34.677431 




55 


C s 


-63.711040 




80 


Cs 


-92.767096 




105 


Ci 


-121.841138 




6 


o h 


-6.873544 


* 


31 


Ci 


-35.838691 




56 


c 2 


-64.871474 




81 


C s 


-93.931045 




106 


C'i 


-123.000416 




7 




-8.028455 


* 


32 


Ci 


-36.999668 


* 


57 


Ci 


-66.031560 


* 


82 


Ci 


-95.092765 




107 


Ci 


-124.164675 


* 


8 


D 2d 


-9.182587 


* 


33 


Ci 


-38.157466 


* 


58 


Ci 


-67.193590 




83 


Cs 


-96.259818 


* 


108 


Ci 


-125.323742 




9 


Civ 


-10.337961 


* 


34 


Ci 


-39.318659 




59 


Ci 


-68.352805 


* 


84 


Ci 


-97.415006 




109 


C'i 


-126.488046 




10 


D 4d 


-11.492834 




35 


Ci 


-40.483344 




60 


Ci 


-69.515486 




85 


Ci 


-98.582433 


* 


110 


Ci 


-127.654226 


* 


11 


c 2 


-12.648250 


* 


36 


D 2h 


-41.646892 


* 


61 


Ci 


-70.678373 


* 


86 


Ci 


-99.742283 




111 


c 2 


-128.824065 


* 


12 


Ci 


-13.804680 


* 


37 


C s 


-42.803277 




62 


Ci 


-71.838429 




87 


Civ 


-100.912912 


* 


112 


C'i 


-129.985671 


* 


13 


c s 


-14.960307 


* 


38 


Did 


-43.967614 




63 


Ci 


-73.002649 


* 


88 


Civ 


-102.077956 


* 


113 


Ci 


-131.150549 


* 


14 


D 6d 


-16.116921 




39 


Civ 


-45.132103 


* 


64 


Ci 


-74.172601 


* 


89 


C s 


-103.236780 


* 


114 


Ci 


-132.311215 


* 


15 


D 3 


-17.276238 




40 


C a 


-46.291870 


* 


65 


Ci 


-75.330653 




90 


Ci 


-104.397182 




115 


Ci 


-133.475490 


* 


16 


C s 


-18.434583 




41 


Civ 


-47.451405 




66 


Ci 


-76.493660 


* 


91 


Civ 


-105.565322 




116 


C'i 


-134.637941 




17 


Civ 


-19.600675 


* 


42 


D 4 


-48.613389 


* 


67 


Ci 


-77.650700 




92 


Civ 


-106.733237 


* 


117 


Ci 


-135.800525 


* 


18 


Civ 


-20.764865 




43 


Ci 


-49.774701 


* 


68 


Ci 


-78.815370 


* 


93 


C s 


-107.892659 


* 


118 


Ci 


-136.963187 




19 


Civ 


-21.923540 




44 


C s 


-50.937992 


* 


69 


Ci 


-79.982620 


* 


94 


C s 


-109.048698 


* 


119 


Ci 


-138.125231 




20 


Civ 


-23.080283 


* 


45 


Civ 


-52.100346 


* 


70 


Ci 


-81.143066 


* 


95 


C s 


-110.214298 




120 


Ci 


-139.285228 




21 


(7i 
^i 


-24 234793 




46 


c 


-53 258607 




71 


1 


-82 302019 




96 


^ 3 


-111 382493 




121 


Co 


-140 450672 




22 


Ci 


-25.390380 




47 


Cs 


-54.421794 


* 


72 


Ci 


-83.462242 




97 


Cs 


-112.543725 


* 


122 


Ci 


-141.614397 




23 


Ci 


-26.552108 
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they are again exceptionally small. The average devia- 
tion from -Eave compared to the average energy per atom 
is 0.67% for Zn and 0.31% for Cd. For comparison, the 
value is 16% for Lennard- Jones clusters, 2.1% for Gupta 
lead clusters^S and 5.6% for aluminium clusters^ in 
the same size range. 

The implications of these small fluctuations are that 
the differences between the more and less stable sizes 
indicated by Fig. [5] are small, making it more difficult 
to observe any such 'magic' numbers. Interestingly, one 
would expect that the properties of fully disordered clus- 
ters would evolve smoothly with size. By this measure, 
these zinc and cadmium clusters seem to be close to this 
limit. 

We will first examine the zinc clusters in detail, and 
then later look at the relatively small structural difference 
between the two systems. A selection of zinc clusters are 
depicted in Fig. [21 that are either particularly stable or 
have some interesting structural feature. Up to N = 10 
the structure of the smallest zinc clusters are typical of 
what one usually finds for clusters modelled by empir- 
ical potentials. However, instead of structures leading 
up to the 13-atom icosahedra, more open structures are 
then preferred. The example shown in Fig. |3| for Zni3 



can be considered as a polytetrahedral fragment of a 19- 
atom double icosahcdron. It is insightful to examine, why 
this structure is lower in energy than the 13-atom icosa- 
hcdron. The icosahedron is in fact 0.335 eV higher in 
energy, and this stems from an unfavourable embedding 
energy. The better pair energy, because of the icosahe- 
dron's greater average coordination number, does little 
to offset this. As the distances between the central atom 
and the vertices of the icosahedron are 5% shorter than 
the distance between adjacent vertices and p(r) increases 
rapidly with decreasing r, the central atom has an ex- 
tremely high value of p, namely 19.568, that is much 
greater than the optimal value. By contrast, the p values 
for all the atoms in the global minimum are relatively 
close to Pxtai- This is achieved by the atoms with low 
coordination number having shorter average distances. 

This example illustrates two important features of the 
potential. First, low-coordinate surface atoms are not 
disfavoured, in contrast to a potential with a strong pair 
component. Second, it is important that interior atoms 
do not have short distances. As we will see, it is of- 
ten much better for an interior atom to have a coordi- 
nation number larger than twelve since then its nearest- 
neighbour distances will be longer than those between its 
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FIG. 2: Energies of the putative global minima of (a) Znjv and (b) Cdjv relative to Save, a four-parameter fit to their energies. 
£?vc = -1.362UV + 0.1062iV 2/3 - 0.00677V 1/3 + 0.09f3. ESI = -1.16847V + 0.0391iV 2/3 - 0.0143iV 1/3 + 0.0425. 



neighbours. This latter feature is illustrated by the 14- 
atom global minimum. The 13-coordinate atom has an 
average nearest-neighbour separation of 1.052 ro, whereas 
the average separation for the nearest-neighbour contacts 
not involving this atom is 0.958 ro. Consequently p for 
the central atom is only 10.381 

The global minima for iV=16-21 can all be charac- 
terized as distorted decahedral structures. Decahedral 
structures are based on pentagonal bipyramids (hence the 
name) and have a single five-fold axis of symmetry. The 
pentagonal bipyramids themselves are usually not partic- 
ularly stable because they have a relatively non-spherical 
shape. By the introduction of reentrant grooves at the 
five equatorial vertices that are parallel to the five-fold 
axis, more stable Marks decahedra can be produced^lS 
An 18-atom example is shown in Fig.0]that was derived 
from a 23-atom pentagonal bipyramid, where the intro- 
duction of the grooves gives rise to five four-coordinate 
capping atoms. 

This 18-atom Marks decahedron is in fact itself not 



that unstable and is only 0.042 eV above the global min- 
imum. Even though the central atom is 12-coordinate, 
this atom is able to maintain a reasonable p value (9.867), 
whilst an expansion of the structure along the axial di- 
rection allows the surface distances to contract. 

However, the energy can be further improved by dis- 
torting this Marks decahedron. The distortion of this 
decahedra is easy to see for the 17-atom global mini- 
mum. This structure can be derived from the Marks 
decahedron (minus one capping one atom) by a diamond- 
square-diamond rearrangement where the contact be- 
tween the atoms in the groove opposite the missing cap 
atom is broken and a contact between the adjacent cap- 
ping atoms is formed (compare the second views of the 
two structures in Fig. |3| and 0J. During this process 
the two capping atoms are also drawn closer to the cen- 
tral atom to form new contacts, thus giving the central 
atom a coordination number of 14. The planes contain- 
ing the two five-fold rings of atoms that were parallel 
in the Marks decahedron are now splayed apart (see the 
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FIG. 4: A selection of reference structures for interpreting the 
Znjv and Cdjv global minima. 

third view of the Z1117 global minimum). It is also no- 
ticeable that the structure in the region of the distortion 
locally resembles the 13-atom icosahedron with the sur- 
face atoms now six coordinated The second view of Znu 
looks down a two-fold axis of this icosahedral-like region. 

The energetic advantage provided through this 
'groove-bridging' distortion is mainly through the low- 
ering of the energies of the two capping atoms that come 
into contact. The increase in their coordination num- 
ber to 6 gives rise to both an increase in p and a more 
favourable pair energy. The energy of the central atom 
is left relatively unaffected, because the increase in co- 
ordination number is compensated by an increase in the 
average nearest-neighbour separation for that atom. 

The other global minima for iV=16-21 can be also un- 
derstood in terms of this distortion. Zni6 simply has one 
less capping atom than Zni7. For iV=18 and 19 two of 
the grooves are bridged, giving rise to a 16-coordinatc 
central atom. The combined effect of the two distortions 
also causes the contact in the groove adjacent to the these 
two distortions to be broken, thus giving rise to a square 
face, which is itself capped in Znig, and a four- fold axis of 
symmetry (see the second view of Znig). Zn2o and Zn2i, 
like Zni7, are only distorted at one groove; the additional 
atoms fill in some of the other grooves taking the struc- 
ture closer towards the 23-atom pentagonal bipyramid. 

Such distorted decahedral structures have 
previously been seen for a number of metal 
potentials, 10 i 13 i 14 i 21 i 36 i 42 i 43 i 44 particularly those that 
have a tendency to disorder. Indeed, Soler et al. have 
identified such structures as particularly important for 
amorphous gold clusters, and have used them as a basis 
to suggest schemes to produce potentially new magic 
number clusters for larger sizes^i 

As the size of the clusters increase, so of course must 
the number of atoms in the interior. Zn25 is typical of 
most of the clusters with two internal atoms. It can be 
considered as two interlocking distorted decahedra. In 
the view shown in Fig. [3] a Znis-like fragment can be 
clearly seen. The exception to these type of structures 
occurs at Zii28- The six- fold symmetric structure is a 
continuation of Zni4 and has four stacked hexagonal rings 



arranged in an antiprismatic fashion. One would expect 
such a structure to have three internal atoms, but instead 
there is a vacancy in the very centre of the structure. 
Again, this is to avoid internal atoms with values of p 
that are too large. 

Zn3o and Zn34 are typical of structures with three and 
four internal atoms. The overall shape reflects the trian- 
gular and tetrahcdral arrangement of the internal atoms, 
but it is again hard to detect any overall order. Zn36 pro- 
vides an interesting exception. The surface atoms have 
basically the same geometry as a recently identified 38- 
atom polytetrahedral structure iiLSH However, there are 
only four (not six) internal atoms, and these are arranged 
as a planar rhombus. 

First at A^=35 and then from N=38 a new series of 
distorted decahedral structures begins. Again they are 
based on a Marks decahedra that has been created from 
a pentagonal bipyramid (in this case with 54 atoms) by 
introducing reentrant grooves that are one layer deep. 
A similar mechanism of distortion is again seen. Atoms 
that are on the equatorial edges of the cluster either side 
of a groove are drawn closer together, although in this 
case a contact is not actually formed — the distortion is a 
diamond-square rather than a diamond-square-diamond 
process. The clusters at iV=35, 39, 41, 45, 47 and 49 in 
Fig-ODprovide examples of structures in this sequence. In 
most of the examples, two such distortions are present at 
adjacent grooves, and so the whole edge between these 
grooves is brought closer to the centre of the cluster, 
forming new contacts with the closest interior atoms. As 
can be seen for the second view of the complete 49-atom 
Marks decahedron, this lead to a splaying of the planes of 
pentagonal rings that would otherwise have been parallel. 

For these larger clusters, these distortions represent 
a smaller perturbation on the overall structure than for 
iV=16-21. Indeed, they are barely visible when the clus- 
ter is viewed down the quasi-fivefold axis. All that can 
be seen is a slight opening of the angle at the reentrant 
groove. 

Comparing the energetics of the distorted and undis- 
torted 49-atom Marks decahedron shows that the main 
improvement arising from the distortion is for the interior 
atoms. The movement of the equatorial edges involved 
in the distortion closer to the centre of the cluster com- 
presses the rest of the surface somewhat and increases 
the coordination number of some of the interior atoms. 
This leads to a structure where the ratio of the nearest 
neighbour distances for the interior atoms to those of the 
surface atoms increases from 1.049 to 1.060. Thus, there 
is less need for the surface to shrink inwards and com- 
press the core in order to improve the p values of the 
surface atoms. Consequently, pbuik, the average p value 
for the interior atoms, decreases from 9.197 to 8.756, i.e. 
closer to p x tai, and so the energy of the interior atoms de- 
creases. By contrast the distortion has little overall effect 
on the energy of the surface atoms. Generally, the atoms 
in the equatorial plane, especially those close to the dis- 
tortion, improve their energies, but many of the other 
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surface atoms lose out due to the breaking or stretching 
of contacts in the axial direction. 

This decahedral series of structures is interrupted at 
N=42 and 43 by structures with four-fold symmetry. 
These clusters are loosely related to the 44-atom fee oc- 
tahedron, but with one or two opposite vertices removed. 
The top half of the cluster is then twisted with respect 
to the bottom to give a slightly buckled outer layer. 

Beyond 7V=50 these groove-bridged Marks decahedra 
are no longer most stable. Instead, there is a size range 
where the structures typically have no overall order, but 
where motifs that resemble fragments of Marks decahe- 
dra and Mackay icosahedra are evident. Zn 54 provides a 
typical example, whereas Zns2 is a structure with more 
order evident. From one side Zn 52 has perfect Mackay 
icosahedral order. What has occurred to the other side 
can be determined by comparing to the 50-atom incom- 
plete Mackay icosahedron shown in Fig. 0] The two con- 
tacts that complete the five-fold rings around the two 
empty vertices of the Mackay icosahedron have been bro- 
ken, removing some of the tension in the outer layer of 
the Mackay icosahedron that makes it so energetically 
unfavourable. 45 Then two atoms (oriented horizontally 
with respect to the views in the figures) are inserted 
into the coordination shell of the central atom, increas- 
ing its coordination number. Such a structure has previ- 
ously been located for AU52 modelled by a Sutton-Chen 
potential^ 

From the particularly stable structures based on the 
18-atom and 49-atom Marks decahedra, it is not surpris- 
ing that there is another series of magic numbers leading 
to the 100-atom Marks decahedron. Again this Marks 
decahedron is formed from a pentagonal bipyramid by 
the introduction of grooves of depth one layer. More gen- 
erally, one can use this trend to predict potential magic 
number clusters outside the size range of this study. The 
sizes for such complete Marks decahedra are given by 

r. gi 

N = -n 3 + 5n 2 + —n + 2 (14) 
6 6 

where n is the number of atoms on an equatorial edge of 
the Marks decahedron. This gives N=1S, 49, 100, 176, 
282, . . . 

The first of this set of groove-bridged decahedral struc- 
tures occurs at A=75 and the last at iV=101. Rep- 
resentative examples at N=75, 77, 79, 83, 86-88, 92, 
96 and 100 are shown in Fig. [3] At the smallest sizes 
the decahedra are still quite asymmetric. From Fig. [21 
one can see that the particularly stable sizes occur at 
N = 100 — 4m, where m=0-3. These structures can be 
formed from Znioo by the sequential removal of 4-atoms 
from each groove to give structures that have m grooves 
that are two layers deep. It is interesting to know that 
the complete Marks decahedron is not in fact the most 
stable of these structures, but that a slight asymmetry is 
preferred. 

Comparing the energetics of the 100-atom distorted 
and undistorted Marks decahedra reveals a similar story 



to Zn4g. The main stabilization of the distorted struc- 
ture is due to the decrease in pbuik from 8.875 to 8.789. 
It is also interesting to analyse the reasons for the sta- 
bility of structures at JV=88 and 92 with four-coordinate 
surface atoms. Firstly, this is because a four-coordinate 
atom can compensate for its low coordination by having 
very short nearest-neighbour distances, thus achieving a 
reasonable p value; e.g. 7.519 for Zng2. Secondly, this 
additional short contact reduces the need for the four 
surface atoms in contact with the adatom to shrink in- 
ward and compress the adjacent interior atoms. For Zng2 
the p values for these atoms only increase by 0.313, de- 
spite this extra contact, because the pair separations for 
the other six contacts have expanded by 2.5%. 

The structures at A=63-70 are also closely based on 
the 100-atom Marks decahedron, but with a further dis- 
tortion. From the first view of Zn^ one can see the 
resemblance to the 63-atom structure in Fig.0]that is an 
asymmetric fragment of the 100-atom distorted Marks 
decahedron. The second view of Zn64 shows the 2-fold 
axis. From the equivalent view of the 63-atom structure, 
it is clear that an extra atom has been added to the col- 
umn of three capping atoms (at the bottom of the first 
view and in a vertical line in the centre of the top half 
of the cluster in the second view), and then a small twist 
has been given to the two halves of the structure. Zngg 
has basically the same structure but with an increase in 
the size of one of the decahedral faces. 

Based on the favourability of small distorted Marks 
decahedra (similarly to those we see at N= 16-21) for 
gold clusters modelled by a Gupta potential, Soler et al. 
suggested a means for using this structural pattern to 
generate potentially stable large clusters^ A 55-atom 
Mackay icosahedron can be generated from a 13-atom 
icosahedron by the addition of atoms at the centre of 
each nearest-neighbour contact. By the same process, an 
80-atom structure can be generated from an 18-atom dis- 
torted Marks decahedron with one of the grooves bridged. 
Soler et al. envisaged that this structure would be most 
stable when the three four-coordinate atoms in this 80- 
atom structure were removed^i However, the suggested 
structure is actually the global minimum for Znso, re- 
flecting the stabilization of low-coordinate atoms that is 
possible for this potential. 

An alternative way to generate the Zngo structure is to 
introduce reentrant grooves two layers deep into the 100- 
atom pentagonal bipyramid. This Marks decahedron is 
then distorted to produce four new contacts that bridge 
the two capping square pyramids adjacent to a groove. 
This distortion is analogous to that for bridging a groove 
one layer deep. In the region of the distortion, the struc- 
ture looks locally like a 55-atom Mackay icosahedron. 

Bridging a groove that is two layers deep is generally 
less favourable than bridging one that is one layer deep. 
There are only five examples for N <100. The others 
occur for Zn7i_73 and Zn§4. Zn7i and Zn§4 are quite in- 
teresting because they involve distortions of grooves that 
are both one and two layers deep. For example, Zng4 can 
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be formed from Zrigo, by first filling in one of the grooves, 
so that it is then only one layer deep and then applying 
the distortion to this groove. 

After the completion of the previous decahedral se- 
ries associated with the Marks decahedra at N=18 and 
49, there are size ranges where the majority of clusters 
have no discernible overall order at 7V=22-34 and N=51- 
62, before structures based on the next Marks decahe- 
dron become lowest in energy. However, even though the 
next complete Marks decahedron is at JV=176, beyond 
AT=101 structures based on this larger Marks decahe- 
dron are immediately lowest in energy. From iV=102- 
117 (excepting iV=108) the zinc clusters are asymmetric 
groove-bridged decahedra with the same twist distortion 
as that for iV=63-70. Like for Zn^ this leads to struc- 
tures with a two- fold axis of symmetry at N=1U and 
115 (see the second views of these structures in Fig. @J. 

The first global minima that are fragments of the 176- 
atom Marks decahedron with only distortions bridging 
grooves one layer deep occur at iV=122 and 123. There 
are also a number of structures, where grooves two layers 
deep are bridged; e.g. AT=108, 120, 121 and 125. The 
structures at 7V=108, 121 and 125 are slightly more com- 
plicated, because the reentrant {111} faces in the grooves 
of the undistorted Marks decahedra are not triangular 
but instead have two atoms along the outer apical edge. 
This is illustrated for Znios, which is based on the Ill- 
atom Marks decahedron shown in Fig. 01 On bridging 
this groove, it is favourable to remove three atoms that 
lie on the ah mirror plane of the Marks decahedron, and 
to twist the top and bottom of the cluster slightly to 
remove a small {100}-type face that would otherwise re- 
sult. 

As has already been noted, for just over half the cad- 



mium clusters in the size range considered, the global 
minimum have the same structures as the zinc clusters. 
A selection of examples where the global minimum differs 
from zinc are depicted in Fig. The global minima are 
again dominated by distorted Marks decahedra, although 
the positions of capping atoms (e.g. Cdi6 and Cdig) or 
the number of grooves that are bridged (e.g. Cdioo) may 
be different. To take Cdioo as an example, this structure 
is second lowest in energy for zinc, because, although it 
has a better embedding energy than the global minimum, 
this is more than offset by a worse pair energy. However, 
because of the reduced magnitude of the pair interactions 
for cadmium this structure is now lowest in energy. 



Another difference between the two systems is the more 
pronounced stability of the twisted Marks decahedra for 
cadmium, as indicated by Fig.|2]and the larger size ranges 
(^=63-72 and 102-121) for which these structures are 
most stable. A further difference is that decahedral struc- 
tures where grooves of depth two are bridged are less 
common — there are only three examples for cadmium 
clusters in the size range considered here. 



The least coincidence between the zinc and cadmium 
global minima probably occurs in the size windows be- 
tween the series of distorted Marks decahedra where the 
structures often have no overall order. This probably 
reflects the large number of disordered structures that 
only have small differences in energy. Examples of clus- 
ters from these size ranges with ^=26, 30, 38, 55 and 56 
are shown in Fig. [5J 
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IV. CONCLUSION 

I have analysed the structures of the global minima up 
to N=125 for two metallic potentials that have been pre- 
viously found to have a particularly strong tendency to 
exhibit disordered clusters. This study confirms that the 
clusters exhibit none of usual ordered forms for materials 
that are close-packed in bulk. Instead, the majority of 
the clusters are based on distorted oblate Marks decahe- 
dra, but where the distortions are well-defined. There are 
series of structures associated with the 18-, 49-, 100- and 
176-atom complete Marks decahedra. However, in be- 
tween there are size windows where the majority of the 
clusters have no overall order, except for N > 100 where 
there is a direct transition between structures based on 
the 100- and 176-atom Marks decahedra. 

Because the effective pair potentials of these zinc and 
cadmium potentials are very shallow, thus providing rel- 
atively little constraint on the pair separations, there is 
a strong tendency for the surface atoms to contract in- 
wards in order for these atoms to obtain a better many- 
body embedding energy. The resulting compression of 
the interior of the cluster causes the conventional struc- 
tures to be disfavoured jil because it results in a Pbuik 
that is significantly larger than p x t a i- Instead structures 
for which the nearest-neighbour distances for the surface 
atoms are naturally longer than those for the interior 
atoms are likely to be favoured, because this allows the 
difference between pbuik and p SUT f to be reduced, bringing 
them both closer to the optimal value. The distortions 
of the oblate Marks decahedra help to achieve just this. 

This paper should be seen in the context of a growing 
research program that has sought to understand how the 
observed structures for metal clusters depend on the form 
of the many-body potentialpii*i 7 *22^£^Zi Even for simple 
potentials of the embedded-atom (or glue) form (Eq. [2J 
this leads to significant complexity. The current results 
further highlight how many-body potentials can lead to 
the stabilization of unusual structural forms iiiL^LSili 
Of course, at sufficiently large sizes the bulk structure, 
in this case hexagonal close-packed, should become most 
stable, but the energy differences found in this study sug- 



gest that the current clusters are far from this limit. 

The emphasis in this paper has been on the zinc and 
cadmium clusters as model systems with a strong ten- 
dency to disorder. As noted in Sect. Ill Al with the Gupta 
potential it is difficult to capture the large c/a ratio for 
Zn and Cd without introducing discrepancies for other 
properties. I have also performed a brief set of global 
optimization runs using Cleri and Rosato's alternative 
Gupta potential for Cd, where the c/a ratio was allowed 
to vary during the fitting. For the selection of sizes I 
tested the usual ordered forms, i.e. close-packed, icosa- 
hedral and decahedral, were always most stable. Reyes- 
Nava et al. also found these cadmium clusters to have 
strong first-order- like melting transitions, again reflect- 
ing the strong ordering for this alternative potential*^ 
Given this strong dependence on the potential parame- 
terization, one should be somewhat sceptical about the 
abilities of these potentials to provide realistic models 
for zinc and cadmium clusters. However, the tendency 
to disorder for these clusters seems to be reproduced by 
electronic structure calculations using density-functional 
theoryi 

There has been little other relevant work. Experiments 
on these clusters have focussed on their electronic shell 
structure42iS£ rather than their geometries, and their are 
only a few theoretical studies to go beyond very small 
sizes. For cadmium clusters electronic structure calcula- 
tions have been performed up to 20 atoms^i These found 
Cdi3 to be an icosahedron and the larger clusters to be 
based upon this structured Ramprasad and Hoagland 
studied zinc clusters using a many-body potential of the 
embedded-atom form (Eq. |2J) , but where the effective in- 
teractions looked very different from the current Gupta 
potential.^ Of the series of candidate structures that 
they reoptimized, they also found icosahedral structures 
to be lowest in energy. 
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